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Abstract The present paper describes the experimen- 
tal implementation of a measuring technique employing 
a slowly moving, near resonant, optical standing wave 
in the context of trapped ions. It is used to measure 
several figures of merit that are important for quantum 
computation in ion traps and which are otherwise not 
easily obtainable. Our technique is shown to offer high 
precision, and also in many cases using a much simpler 
setup than what is normally used. We demonstrate here 
measurements of i) the distance between two crystalline 
ions, ii) the Lamb-Dickc parameter, iii) temperature of 
the ion crystal, and iv) the interferomctric stability of 
a Raman setup. The exact distance between two ions, 
in units of standing wave periods, is very important for 
motional entangling gates, and our method offers a prac- 
tical way of calibrating this distance in the typical lab 
situation. 
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1 Introduction 

During the past decade ion traps have been shown to be 
one of the most promising candidates for the implemen- 
tation of quantum computers [I] and quantum simula- 
tors PJ. A basic building block for these two scientific 
endeavours is the controlled interaction of several ions 
with an optical standing wave. In the case of quantum 
computing, a tool for creating the necessary entangled 
states is the controlled motion of the ions that will couple 
via the Coulomb interaction, such as in the case of the 
geometric phase gate j3]. In this type of gate, a specific 
collective vibrational mode is chosen to mediate the cou- 
pling, high fidelity interactions is possible only with pre- 
cise control of the experimental situation. This includes 
good knowledge of parameters such as the alignment of 
the ions in the standing wave and their localization (ion 
crystal temperature) , in order to reside within the Lamb- 
Dicke regime where the ions are localized better than the 
wavelength of the standing wave. These parameters can 



be very challenging to actually measure. Fig. [T] depicts 
a typical experimental situation with an ion crystal in a 
standing wave from two laser beams irradiated from out- 
side the apparatus. The actual wavelength of the beat 
pattern along the trap axis depends on the angles of the 
beams with respect to this axis and is therefore hard to 
set a priori. As the ion distance in a typical experiment 
is 10 1 — 10 2 times the beat pattern period, an alignment 
precision of better than 1% is clearly required. Here we 
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Fig. 1 Alignment of the two components of the beams to- 
gether with the resulting direction of the Ak-vector (given 
by km — ki^) with respect to the trap axis. 



present a new technique, which is based on a slowly mov- 
ing standing standing wave, designed to measure several 
important parameters for controlling entangling opera- 
tions in ion traps. We show measurements of i) the dis- 
tance between two crystalline ions, ii) the Lamb-Dickc 
parameter, iii) the temperature of the ion crystal, and 
iv) the interferometric stability of the Raman laser beam 
setup that is used to excited the ion motion. The basic 
idea of the standing wave technique is to use two laser 
beams, tuned close to a resonance of a dipole transition 
that can both cool the ions and yield fluorescence, but 
with a very small frequency offset with respect to each 
other, in our case chosen to be 2 Hz. This creates a very 
slowly moving wave, in which ions in a peak of the wave 
yield fluorescence but ions in a valley does not. For a 
two ion crystal in particular, a strong 2 Hz fluorescence 
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beat pattern will be detected only if both ions are in 
peaks of the wave simultaneously, and thus the ampli- 
tude of the beat pattern becomes a direct measurement 
of the distance between the two ions, in units of the 
wave period. The technique is shown to offer very good 
accuracy, allowing measurements of the Lamb-Dicke pa- 
rameter in the sub-percent range. Looking at this from 
the opposite point of view, one can also say that the 
ions are used to map the properties of the light with a 
sub-wavelength precision, similarly as to what has been 
discussed before [1] . From this perspective, we are oper- 
ating in a range between the single ion resolution of [3] 
and the completely averaged regime that is obtained for 
large ion crystals [S|. The method also operates with- 
out ground state cooling, which together with the fact 
that the beating pattern that can be observed in real 
time, enables the method to be used as a practical ion 
distance calibration tool. In the context of probing con- 
densed matter objects by light or particle scattering, our 
method represents a measurement of the structure factor 
and Debyc- Waller factor of the ion crystal, although with 
a fixed wavenumber of the probing light and a structure 
that can be varied by means of external control param- 
eters. This could be used as an ideal tool for studying 
implementations of the Frenkel-Kontorova model with 
trapped ions as proposed in Ref. [BJ. 

The paper is organized as follows. In Sec.O the tech- 
nique is explained in detail, with the experimental setup 
described in Sec. [31 and with the results of measuring 
the different parameters presented in Sec. |4] In Sec. [3 
we elucidate in detail why and how this technique can 
be useful for the optimization of entangling gate opera- 
tions, finally in Sees. [BJ and [7J we discuss the suitability 
of our method as a standard measurement tool and a 
possible generalizations to more ions. 



2 Methods: Standing wave technique 

For a single ion, the fluorescence intensity from two laser 
beams with the same polarization and with a detuning 
of Alo with respect to each other is given by 



wq is given by 



h{t) = A + B cos (Ak cS ■ x - Aujt + A<j>) 



(1) 



where Ak e g is projection of the difference k-vector be- 
tween the two beams onto the trap axis and A<p is their 
optical phase difference. The term A is the background 
fluorescence signal, which will be ignored in the follow- 
ing. The amplitude of the temporally varying signal B 
depends on the laser power and the specific saturation 
of the optical transition as well as on the overall detun- 
ing of the two beams relative to the line width. The laser 
power is chosen such that the transition is not saturated. 

The distance of two crystalline ions in a common 
harmonic potential characterized by the trap frequency 
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For two ions with a center of mass position xq, we can 
describe the beating component of the fluorescence as 
sum of the fluorescence rates from each ion: 



B cos(Z\fc cff • (x + l ) - Aut + A4>) + 
B cos(zifc c ff • (xq — Iq) — Aut + A(j>), 



(3) 



which is implicitly dependent on Xq and lo. Note that 
this signal is not due to interference of fluorescence from 
the two ions [7J, which other imaging schemes use |S|, 
but rather due to the fact that the fluorescence from 
the two ions is correlated or anticorrelated depending in 
their alignment with respect to the standing wave and 
that the detection does not provide spatial resolution. 
In the case of a detection apparatus providing spatially 
resolved detection of ions [3], cross correlation signals 
can be used in the same manner. 

We need to take into account that the ions are delo- 
calizcd due to thermal and quantum fluctuations of their 
position. The probability of finding the ion at a ccntcr- 
of-mass position x is given by a Gaussian probability 
distribution^: 
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The width of the thermal distribution, ctoi depends on 
the mean phonon number 



n 



(6) 



(o) 



and on the width of the ground state wavefunction, a 

2mwo • From this it is clear that the width of the 
thermal packet depends in two ways on the axial trap 
frequency luq: firstly because a change to the trap fre- 
quency signifies a change to the confinement which also 
changes the width of the ground state packet, and sec- 
ondly because a different trap frequency also gives a dif- 
ferent mean phonon number for a fixed amount of ther- 
mal energy. A similar argumentation holds for the fluc- 
tuations of the ion distance Iq + I associated with the 
stretch mode vibrating at oj\ = \fiuiQ. We now integrate 
over all fluctuations, such that the final fluorescence in- 
tensity of the beating component as a function of time 
and ion separation is given by 



h(t) = II I 2 (t)po(x') Pl (l')dx'dl' 



(7) 



1 Note that a x = ao/y/2 



Precision measurements in ion traps using slowly moving standing waves 



3 




Time (s) Frequency (Hz) 

Fig. 2 (color online) signal beating taken at two different ion distances, controlled by voltage on the segments of the microtrap. 
Part a) and c) show the raw data (blue curve) taken at a peak and a bottom of the oscillation strength respectively. Part b) 
and d) show the corresponding Fourier power spectra. The amplitude of the 2 Hz Fourier component has also been included in 
the left part of the figure in red. The strong oscillations correspond to the point at 4.82 [im in Fig. fj] and the weak oscillations 
to the point at 4.73 \im. 



where I2 (t) is the two- ion intensity from Eq. using the 
fluctuating coordinates with their corresponding proba- 
bility distributions from Eq. ((4]). This integral yields 

h(t) oc cos(A^) cos(Ak cS l )e- {1/4)Ak ^( a ° +a *) , (8) 

where x$ = and A<p = A(f>(t) represents phase insta- 
bilities between the interferomctric branches. This ex- 
pression shows that we obtain a signal on the beating 
frequency, Auj, which has an amplitude that depends on 
the ion separation. As can be seen, when the ion sep- 
aration is exactly an integer number of wavelengths of 
the k-vector, the beat amplitude has a maximum, and 
for half-integer numbers, the beating disappears entirely. 
The beat amplitude also depends on the temperature of 
the ion crystal in the form of the wave packet width from 
the two different vibrational modes, gq and G\. 

3 Experimental setup 

The ions used in these experiments were 40 Ca + , trapped 
in a segmented micro-structured Paul trap [10] , with ra- 
dial trap frequencies of of cj la( j/27r = {2.4,3.0} MHz. 
The axial dc trapping voltage is varied during the ex- 
periment in order to change the distance between the 
ions. This causes the axial trap frequency u>o to vary be- 
tween 2-7T ■ 0.96 and 27r • 1.4 MHz. The main laser inter- 
action during the experiment consists of a pair of beams 
at 397 nm, both linearly polarized, hitting the ion at 
approximately 90° angle with respect to each other, as 



indicated by Fig. [TJ These beams have a small detun- 
ing from each other of Aoj = 2ir ■ 2 Hz, controlled by a 
waveform generator. They are also red-detuned slightly 
from the optical transition, such that they can provide 
the necessary Doppler cooling to the ions and serve as 
the source for the fluorescence that is detected. In ad- 
dition, there are also lasers on 866 and 854 nm for the 
rcpumping purposes. The beating fluorescence light can 
be monitored in real time by a CCD camera, but for 
technical reasons, a photo multiplier tube was used to 
capture the data presented in Sec. [U A binning time 
of 50 ms was used, and the measurement durations for 
each run was 50 s. There is also an external magnetic 
field aligned at 45° with respect to the trap axis for the 
purpose of breaking the degeneracy of the ground level 
spin states. 

4 Results 

Measurements of the fluorescence were taken at a total 
of 33 different ion distances, with 5 runs of 50 seconds 
detection time for each distance. Examples of such a 
runs, when the ion distance was corresponding to a peak 
and a valley of the beating respectively, are displayed in 
Fig. [5J For proper analysis of the beating, we perform a 
Fourier transform of the signal. The Fourier spectrum is 
then multiplied with a supcrgaussian window 

= e - < ^" ,2 \ (9) 
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where the center frequency is chosen to only keep the 
beating component, Auj = 2n ■ 2 Hz, and where the order 
was chosen to be n = 4, giving a window that is roughly 
flat around the interesting parts, and zero elsewhere. The 
FWHM width of the window was chosen to be 0.3 Hz, 
but other widths was tested and the result was found to 
be relatively independent of the choice of window width. 
If one multiplies the window function to the positive 
side of the Fourier plane only, then upon inverse Fourier 
transform, a complex function, A{t)e~ l ^ lyt \ is obtained. 
Here, A(t) corresponds to the amplitude of the 2 Hz 
component of the signal as a function of time, and is 
plotted as the red ovcrlayed curve in Fig. [2] Similarly, the 
phase factor, 4>(t) corresponds to the phase evolution of 
the 2 Hz component and it plotted in Fig. [3J This figure 
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Fig. 3 (color online) phase angle of the 2 Hz beating com- 
ponent as a function of time is plotted as the black (solid) 
line with corresponds axis to the left. Its derivative is plot- 
ted in red (dashed), with axis to the right, which gives a 
measurement of the frequency stability of the interferometric 
setup. 



also contains the derivative of the time dependent phase 
angle, which gives a measure of the frequency stability 
between the two arms of the interferometric setup. From 
an average of the standard deviation of the derivative 
for several traces we find an interferometric coherence 
time of about r ~ 12 s. The interferometric stability is 
important for other types of experiments, as this pair of 
laser beams is the same as is used to performed Raman 
type quantum gates in the quantum computer scheme, 
and we see that the coherence time is much longer than 
the typical gate times, which are in the order of 10 to 50 
us. We therefore expect that phase fluctuations during 
such gate operations are much less than 7r/100, thus not 
limiting the gate fidelity. Interferometer stability might 
also be important when applying optical light forces with 
fixed phase relations to other forces, such as electric drive 
forces from the trap electrodes, though this was not used 
in the context of this paper. 

Since the Fourier operation has picked out only the 
amplitude part of the beating, we can now write Eq. ([5)1 
as 

h{t) cx cos(Aji) cos(Z\fc cff Zo)e- (1/4)4fcrff (" o2+ ^) . (10) 
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Fig. 4 (color online) shows the measured beat power of the 
2 Hz component of the fluorescence as a function of the ion 
distance, together with a theoretical fit using Eq. (|10jl . From 
the oscillation period, the effective wavelength of the beat 
pattern can be determined to A = 267.8(2) nm. Similarly, the 
decay of the oscillations for larger distances (lower trap fre- 
quencies) gives an ion crystal temperature of about 3.7) mK. 



In Fig. |4j the intensity of the beating component has 
been plotted as a function of ion distance for 33 dif- 
ference distance. Eq. (fTTJj) has then been used to make 
a fitting to the result using the overall amplitude scal- 
ing together with Ak and the ion crystal temperature, 
which gives both a x and <j;, as fitting parameters. De- 
spite the sizable vertical error bars for each data point, 
the wave period can be determined with sub-percent ac- 
curacy, yielding a value of the standing wave wavelength 
of A = 2n/Ak cS = 267.8(2) nm. Note that the outstand- 
ing accuracy comes from the fact that there are many 
wave periods between the ions, and as such, even a small 
change to the wavenumber will be multiplied and change 
the fit substantially. The temperature was found to be 
T = 3.7(2) mK, where the accuracy is in the 5% range. 
Using this temperature, the common and stretch mode 
vibrational quantum numbers can be determined from 
Eq. ©, and we find no ~ 56(4) and n\ sa 32(2), for a 
trap frequency of ujq = 1.24 MHz, which corresponds to 
the same ion distance that was chosen for Figs. [2] and [3] 
These values are consistent with the expected results 
of Dopplcr cooling and with measurements using other 
techniques |11| . Furthermore, for this trap frequency and 
with the value of k above, we calculate the Lamb-Dicke 

parameter to be r/ = k c g 



0.237. 
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5 Relevance for High Fidelity Entangling Gates 

The main motivation for the determination of the ion 
distance in units of the standing wave period is the re- 
quirement for high-fidelity quantum gates. The most com- 
monly used entangling gates are based on the application 
of forces mediated by a spatially and temporally varying 
laser beat pattern, where the sign of the force can de- 
pend on the joint internal state of the ions. Two widely 
used gate schemes are the geometric phase gate [5] and 
the Molmer-Sorensen gate P^llTBIHl] . In the first scheme, 
the ion motion is near-resonantly driven by an oscillating 
ac- Stark shift arising from an off-resonant laser field; in 
the second scheme, simultaneous driving of red and blue 
sideband transitions gives rise to an effective entangling 
spin-spin interaction. In both cases, the driven collective 
vibrational mode is displaced and restored to the origin 
after the gate time, leading to the pickup of a geometric 
phase which is conditional on the joint internal state of 
the ions. The relevance of a precise alignment of the ions 
in the driving field is illustrated in the following for the 
case of the geometric phase gate driven by a bichromatic 
off-resonant laser field: The classical force on one specific 
ion, i, is given by its position, x, in the standing wave: 

Fi(x t ) = Asm(Ak e «Xi - St + A4>)m { ^ ', (11) 

where m^) is the spin quantum number for ion i [mj 
for species without nuclear spin), Ak e s is the effective 
wavevector along the trap axis, and A<f> is the phase dif- 
ference of the two beams. Aui is the difference between 
the two-photon detuning and the vibrational frequency, 
i.e. the effective detuning of the drive from the collective 
vibrational frequency. The total force driving the STR 
mode of a two ion crystal is the given by the difference 
of the individual forces acting on the ion, which is due 
to the fact that the forces on the ion need to be opposite 
for driving the STR mode. The spin configuration of the 
crystal in conjunction with the ion distance now deter- 
mines the acting drive force, e.g. if the laser beams are 
tuned to drive the STR mode, only the even spin configu- 
rations (|tt) an( A 144-)) are driven by the laser wave if the 
ion spacing (Eq. [2]) corresponds exactly to a half integer 
multiple of drive period, 21q = (n+l/2)A e g. By contrast, 
only the odd components are displaced if the spacing is 
exactly an integer multiple of the drive period. For a con- 
ditional phase gate, the drive amplitude such that the 
desired differential geometric phase between even and 
odd state is picked up, A<1> = <P cvcn — <P dd- In the case 
of nonidcal ion spacing, corresponding drive amplitudes 
and durations can still be found if the difference of the 
forces acting on odd and even states is large enough. 
However, nonidcal ion spacing will lead to a larger re- 
quired pulse area. The minimum possible amount of laser 
power can therefore only be used upon precise calibra- 
tion of the ion spacing. Any increased amount of laser 
power will lead to increased decoherence by phase and 



intensity fluctuations. For off-resonant gate interactions, 
stimulated Raman scattering will lead to a decoherence 
proportional to the total employed pulse area [15l[T6j 
even in the limit of technical perfection. The precise dis- 
tance calibration becomes even more important for fast 
gates, when stronger drive amplitudes are used and the 
displacements become large enough such that effects be- 
yond the Lamb-Dicke regime have to be considered |TTJ 
[TT] . Here, the precise knowledge of the driving strengths 
is important to chose the gate time such that the dis- 
placement is completely undone at the return time and 
no information about the spin is left in the motional de- 
gree of freedom. Envisaged gate schemes for larger ion 
numbers enabling quantum computing in decoherence 
free subspaces have more stringent requirements on the 
crystal alignment |18j , which will be covered in the next 
section. Furthermore, closed loop strategies such as op- 
timal control theory might turn out to be helpful for the 
realization of quantum gates |19j . however, they require 
the system parameters as input such that the knowledge 
of the effective wavelength is would even be of crucial 
importance. 

6 Generalization to an arbitrary ion number 

We now consider the magnitude of the fluorescence os- 
cillation signal for a linear ion crystal consisting of an 
arbitrary number, N, of ions. The position of the z-th 
ion is given by 

N 

Xi = xf> + ^ a ijQj ( 12 ) 
i=i 

with the equilibrium positions xf^ , the normal coordi- 
nates qj pertaining to the N axial normal modes and the 
Hessian matrix Ref . [5U] . The oscillatory compo- 

nent of the fluorescence signal for ion i is given by 

/W = COs(Aujt) COS (Ak e ffXi) 

= Re (e iAk *« x * 0) e lAKn ^ aim ) , (13) 

where the reflection symmetry of the crystal about x = 
has been used and where Acf> = as well as a homoge- 
neous light field intensity has been assumed without loss 
of generality. Averaging over the thermal and quantum 
fluctuations yields 

/« {Aoj) = cos(A,i)Re |e^«*f TJ j (14) 
where 

fij = / d 1j e 1 

J —OQ 

= e -(i/2)4fe V^I. (15) 
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Fig. 5 (color online) Theoretical beat signal amplitudes ver- 
sus trap frequency for linear ion crystals of different size. Note 
the periodic structure in the case of four ions, which can still 
be placed such that all spacings are integer multiples of the 
beat pattern period. This is not the case for larger crystals, 
as can be seen for the cases of six and twelve ions, where the 
irregularity increases towards larger ion numbers. 



The total beat amplitude is the given by the sum over 
the ion's individual contributions I n (t) — J2i In (*)• The 
total beat amplitude per ion versus trap frequency is 
plotted in Fig. [5] for different ion numbers. It shows 
that the individual ion localization is better for larger 
ion numbers. This counterintuitive behavior can be ex- 
plained by the fact that for the low-frequency eigen- 
modes, the contributions the each ion's delocalization 
o>ij are smaller. In case of the COM mode, this is sim- 
ply due to the fact that the effective oscillator mass is 
n-times the mass of each ion. The high frequency modes 
are subject to a lower initial population according to Eq. 
([5]) , such that they do not contribute significantly to the 
delocalization at sufficiently low temperatures. 



7 Conclusion and outlook 

We have presented a technique where the fluorescence 
from a two-ion crystal in a slowly moving standing wave 
was used to accurately measure the distance between 
ions in units of the standing wave wavelength. Being able 
to set the ion distance exactly to an integer number of 
the wavelength is important for maximizing the fidelity 
of quantum gate operations that are based on exciting 



motional interactions. In addition we show that this flu- 
orescence signal also can be used to obtain the temper- 
ature of the ion crystal as well as information about the 
interferometric stability between the beams. Although a 
figure for the localization of ions in the standing wave, 
i.e. the ion temperature, is obtained, the method is not 
directly suited as a thermometry tool, as the tempera- 
ture is determined in this case by the Doppler cooling 
exerted by the standing wave. The accuracy is limited to 
about 5%, however the method can still be used to assert 
that the localization of the ions, i.e. the quality of the 
Doppler cooling, is slightly above what is expected from 
the Doppler cooling limit. Explanations for this might be 
competing Doppler heating effects on the micro-motion 
sidebands, or uneven cooling depending on sidebands 
and ion position (for further discussion on this, see |21j). 
Since the technique is based on continuous laser interac- 
tion on a slow time scale, monitoring the signal can be 
done in real time without additional experimental obsta- 
cles like ground state cooling or driving coherent dynam- 
ics on motional sidebands. It therefore allows for direct 
alignment of ion crystals to the standing wave by fine 
adjustment of the trap frequency. In the future, we plan 
on applying this method for optimizing the properties of 
four-ion crystals prior to gate operations in decoherence- 
free subspaces [Tg] , 
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